Thermalization of Sputtered Particles as the 
Source of Diffuse Radiation from High 
Altitude Meteors 

Dejan Vinkovic 

School of Natural Sciences, Institute for Advanced Study, Einstein Drive, 

Princeton, NJ 08540, USA 



Abstract 

High altitude meteors become luminous at altitudes above ~130 km, where the 
standard ablation theory of meteor light production is not applicable. The physical 
mechanism responsible for their glow has not been known. We present a model 
that explains their existence, morphology and lightcurves. The model is based on 
particles ejected from the meteoroid surface through the sputtering process. The 
kinetic energy of such a sputtered particle is typically more than 1,000 times larger 
than the energy of particles in the surrounding atmosphere. Thus the sputtered 
particle creates a cascade of collisions in the atmosphere during thermalization. We 
show analytically that this process is capable of producing enough light for detection. 
We also explain the observed relationship between the beginning height of high 
altitude meteors and their maximum brightness. In addition, meteors are modeled 
with a Monte Carlo code developed specifically for this phenomenon. Theoretical 
images reproduce the observed shapes and sizes of high altitude meteors. Their 
exact shape depends on the relative angle between the meteor path and observer's 
line of sight. 
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1 Introduction 

A meteor is a luminous phenomenon produced by a meteoroid colliding with 
the earth's atmosphere. The collision starts with atmospheric particles imping- 
ing on the meteoroid's body, by which the surface heats up and sublimates. 
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The cloud of vapor particles surrounding the meteoroid collides further with 
the atmosphere and becomes ionized and excited. The meteor light is produced 
by deexcitations, with metals from the meteoroid dominating the spectrum. 
When the meteor reaches altitudes where the Knudsen number Kn (ratio be- 
tween the atmospheric mean free path and the meteoroid size) is less than 
~100, vapors completely engulf the body and shield it form direct collisions 
with the atmosphere. After that the heating and evaporation of the meteoroid 
is caused by hot vapor and radiation from the shock front. During the whole 
process the meteoroid body is ablating, that is losing mass either as gas or as 
fragments. The emission of light exists only as a consequence of the ablation 
process. 



The limiting altitude at which most meteo rs become luminous is ~120 km or 
less, with very few startin g at ~ 130 km fe.g. lBetlem et alll200(l l ICampbell et al.1 



20001 : iBrown et all [2002). The atmospheric mean free path is on the order of 
meters at these altitudes, which is enough to induce intense evaporation of 
volatiles from large (centimeters in size) fast meteoroids. It came as surprise, 
therefore, when meteors of much high er beginning heights were detected in 
1995 and 1996 by Fuiiwara et al.l (Jl998) using TV cameras. They reconstructed 
trajectories of two Leonid fireballs of -4 m and -7 m visual magnitude and found 
that they started at ~160 km altitude. This is abo ut 30 km higher tha n the 
highest meteors detected by photographic systems ( Betlem et al. . 2000h . The 
detection was attributed to the much higher sensitivity of TV cameras, espe- 
cially the infrared region of 900 nm. 



This re sult was confirmed and extended further in 1998 by ISpurnv et al 
(|2000a[ ). using high-sensitivity TV systems with limiting magnitude +6.5 m . 
They detected thirteen Leonid fireballs at altitudes above 130km. Three of 
them were above 180km, with the highest at 200 km. The atmospheric mean 
free path is over 200 m at an altitude of 200 km, and more than 50 m at 160 
km, while the meteoroid body has a size of no more than tens of centimeters 
for the brightest detected meteors. It is puzzling how meteors can produce 
light in such conditions. The meteoroid surface is still too cold for sublimation 
or evaporation and, thus, these high altitude meteors can not be explained by 
the standard theory of meteor ablation. 



Spurnv et al. (|2000bh described the morphological features of these meteors. 
Initially they appear as a several kilometers big diffuse spot that rapidly trans- 
forms into a diffuse comet-like V-shape structure with a head and tail. Some- 
times structures like arcs, jets and streamers are visible in the tail. At altitudes 
of ~130 km the meteor changes its appearance into an ordinary droplet-like 
head shape with trail. This also marks the moment when the standard ablation 
theory becomes applicable. High altitude meteors up to 150 km were also de- 



tecte d in other meteor streams: 77-Aquarids, Lyrids, and Perseids ([Koten et al. 
2001k 
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The decreasing size of the diffuse structure with altitude suggests that the 
atmospheric mean fre e path is important as the light is induced by collisions. 
Spurnv et al.l ([2000al ) noted that the beginning meteor height grows with the 
peak brightness of meteors. This implies that the number of collisions increases 
with the meteoroid size. Considering these properties of high altitude meteors, 
we suggest that they are produced by a cascade of collisions in the atmosphere 
during thermalization of particles sputtered from the meteoroid surface. Sput- 
tering is a process in which an energetic particle projectile ejects atoms from 
the surface of a solid material. It has been used in meteor phy sics only to 
study heating of t he surface or slowing down of microme teoroids f Opik . 19581: 
Lebedinets et al. . 19731 : Bronshten . 19831 : Coulsonl . 2002j b The luminosity has 
been associated exclusively with the phase of intense evaporati on 1 , during 
which the number of released atoms is increased dramatically ((Bronshten . 
1983h . Unlike vaporization, sputtering can exist at high altitudes because it 
does not depend much on the meteoroid surface temperature. 



It is the purpose of this work to reproduce basic properties of high altitude 
meteors. We present our 3D Monte Carlo model that traces multiple particle 
collisions within the atmosphere and show modeling results in the form of 
theoretical images from a ground-based CCD detector (TV camera). 



2 Sputtering from the meteoroid surface 



Meteoroids in meteor showers enter the earth's atmosphere with geocentric 
velocity ~15 km/s to 71 km/s. From the point of view of the meteoroid sur- 
face the atmospheric molecules behave like projectiles with kinetic energies of 
hundreds of eV. For example, an oxygen molecule at 71 km/s has 840 eV of 
kinetic energy. Collision with the surface dumps this energy into displacements 
of atoms in the solid until the projectile stops. If the energy is high enough, 
some atoms are ejected from the surface, and we call these sputtered particles. 
If the projectile exits the solid then we call it a scattered particle. 

The sputtering yield Y is defined as the number of sputtered particles rela- 
tive to the number of projectiles. The yield depends on the projectile energy, 
physical and chemical properties of the solid target, and on the type of the 
sputtered atom. The experimental data on sputtering at the energies of inter- 
est here are scarce, thus we mostly rely on numerical models. Since the high 
altitude meteors detected so far belong to cometary material, we can assume 
that they are made of a fragile (low binding energy) amorphous material. 

1 During the preparation of this manuscript, we became aware of a study by Hill, 
Rogers & Hawkes (accepted for publication in Earth. Moon, and Planets) describing 
sputtering as a possible source of luminosity in high altitude meteors. 
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Calculations show that sputtering yields on such targets approach Y=l for 
projectile energies of hundreds of e V and that t he sp uttered atoms have typ- 



ically 5-10% of the impact energy (jField et all 119971 ). In the case of silicates 



typically present in interstellar grain s, the yield can be one or two orders of 
magnitude smaller ( Mav et al. . 2000[ ). Since in this work we do not investi- 



gate details of the meteor spectrum, we do not differentiate between different 
species of sputtered particles. We focus only on the number of particles origi- 
nating from the surface, and therefore describe both sputtered and scattered 
particles with Y. 

The velocity of sputtered and scattered particles has one component equal to 
the meteoroid velocity. An additional (smaller) component is present because 
a particle picks up some additional energy from the impact. The total velocity 
is a vector sum of those two components. This makes it about 100 times faster 
than the average thermal speed in the atmosphere. Such a fast particle will 
create a cascade of collisions before slowing down to the thermal level. Dur- 
ing this process of thermalization, collisions will excite atoms and dissociate 
molecules, which results in production of light by deexcitations. 

We make an estimate of whether this process creates enough light to be de- 
tected by the TV cameras used in the observations of high altitude Leonid 
meteors. During the time exposure At of one video frame, a meteor of ve- 
locity v m travels a distance t> m At. The number of collisions with atmospheric 
particles is then R^iTVmAt n atm , where R m is the meteoroid radius and n atm is 
the atmospheric number density. Each particle originating from the meteoroid 
surface will result in many photons from excited (mostly atmospheric) atoms 
in the cascade. The total number of emitted photons is 

photon = lYR 2 m 7rv m At n atm , (1) 



where 7 is the total number of emitted photons in one cascade. 

These photons are emitted within the overall volume of ~ 47r(L/2) 3 /3, where 
high altitude meteors have size, L, of several kilometers. On the other hand, 
one pixel in the CCD of a TV camera covers only a fraction of the meteor 



angular size on the sky. ISpurnv et al.l (|2000al ) used an imaging setup with 



angular pixel size ~2.5" x2.5". This covers an area, A, of 150x150 m 2 at 200 
km distance. One pixel collects photons from the volume AL. If photons are 
emited uniformly within the volume of the meteor, and asumming a meteor 
size of L=10 km, the fraction of photons contributing to the flux measured 
in one pixel is ^ = 24t4/47tL 2 ~0.0004. However, photons are not emitted 
uniformly, but with a radial gradient. The central region, which we recognize 
as the meteor head, can emit ~100 times more photons than the average, as 
we will see in §4 from numerical models. This puts the estimate on \1/ closer 
to 4%. 
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If the average photon energy is Eq then the measured flux in one pixel is 
p *E Nl& ton 



AnD 2 At 



(2) 



where D is distance to the meteor. To work in units of stellar magnitudes, 
we compare this flux with the flux of Sirius, -Fsirius- Although the exact stellar 
magnitude depends on the spectral sensitivity of the detector, we assume here 
a visual magnitude of —l m for Sirius. The apparent meteor magnitude is then 

m = -1- - 2.5 log f 1I3^E^) (3) 



Spurnv et al. ( 2000a ) and Fuiiwara et al. ( 1998f ) argue that the sensitivity of 
TV cameras to wavelengths of ~800 nm is one of the reasons why they can 
detect high altitude meteors. Photons of that wavelength have Eq=2.5x 10 -19 
J. One collisional cascade can produce 7 ~100 photons (see §4), as a product 
of excitatio n by collisions or diss ociation of molecules. The highest Leonid 
detected bv ISpurnv et al.l (J2000aJ) had _R m =0.07 m and was first detected at 



distance D ~200 km. By using v m =71 km/s for Leonids, n atm ~10 m 
and -Fsirius ~2xl0 -8 W/m 2 (integrated flux within the sensitivity range of 
the camera), our estimate of the brightness is 9Jt ~6.5 m . Thi s is exa ctly the 



limiting magnitude of the video system used bv ISpurnv et al.l (|2000aJ). 



For a given observational setup, the limiting magnitude is, by definition, the 
same as the meteor brightness at the beginning height, H begin . Equation 3 
shows that H begin depends on the meteoroid size. Larger meteoroids reach 
higher maximum brightness <)Jl max during their disintegration below ~90 km, 
thus a correlation between these two parameters is to be expected. Indeed, 
observations show that H begin is higher for brighter meteors (figure 1). In 
addition to the meteoroid size R m and distance to the beginning point D, this 
dependence comes from variables in equation 3 that differ in two meteors from 
the same meteor stream. The size of a meteor scales with the atmospheric mean 
free path Z atm , thus $ oc l~ tm . For altitudes between 160 km and 200 km, the 
atmosphere has / atm oc exp(H begin / 27 .3km) and n atm oc exp(— H begin /26.8km) 
(the U.S. Standard Atmosphere 1976). 

In the single body approximation of ablation, the maximum meteor brightness 
fffimax j n magnitudes is (see Appendix) 



9JT 



—7.5 \ogR m — 2.5 log(cos Z) + constant, 



(4) 



where Z is the zenith angle of the meteoroid's trajectory. Combining this with 
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equation 3 yields the dependence of Hb egin on ^Jl max : 

H begin = -5.55km [Tl max ] + constant, (5) 

where [9Jl maa; ] is the reduced maximum brightness: 

[M max ] = m max + 7.5 log ( ° ) + 2.5 logfcos Z) (6) 

V 200km/ 

scaled to the distance of 200 km. 

Observed values of Hb egin as a function of [9Jt ma:E ] are plotted in figure 2. Linear 
regression gives a slope of -5.9±0.5 km, which agrees with the theoretical slope 
in equation 5. 



3 Monte Carlo model 



Unusual morphology is another characteristic of high altitude meteors. In or- 
der to explain changes in their appearance with altitude, we preformed Direct 
Simulation Monte Carlo (DSMC) modeling. The simulation follows a sput- 
tered particle and all high velocity particles produced by collisions during the 
thermalization process. The three-dimensional distribution of collisions pro- 
vides information on the spatial distribution of photons produced. 

Although DSMC is widely used for modeling interactions betw een rarefied 
gases and solids, including reentry o f spacecraft ( Oran et al. . 1998| ). it has very 



rarely been used in meteor physics. Bovdl ( 2000[ ) modeled the flow around a 1 



cm Leonid meteor at 95 km altitude. The model shows that meteoroid vapors 
surround the body and interact with the atmosphere on a larger scale than 
when the vapor s are not p r oduce d. This confirms theoretical expectations of 



vapor shielding. IZinn et al.l ([20041 ) also performed, but not with Monte Carlo, a 
numerical study of Leonids at these and lower altitudes with a hydrodynamics 
code coupled with radiative transfer. They found that the meteoroid vapors 
are stopped very quickly by the air and moved into the meteor wake, where the 
vapor-air mixture expands and cools. This stopping process is the dominant 
source of energy deposition in the atmosphere. These models demonstrate the 
importance of correct numerical treatment of the vapor cloud. 

All these numerical models have been applied to conditions in which Kn<10, 
typically for altitudes below ~100 km. This is where the intense vaporization 
and compressed air in front of the meteoroid create conditions appropriate for 
the approximation of vapor shielding and the continuous flow regime. This 
choice of flow conditions is also influenced by the numerical codes used in 
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simulations. They are not designed specifically for meteors and their applica- 
bility requires the hydrodynamic fluid regime. In contrast, our DSMC code is 
written and designed specifically for meteors under the molecular flow regime 
at higher altitudes. 

In this study we focus solely on the sputtering effect at altitudes above ~130 
km as relevant to high altitude meteors. For simplicity, we do not differenti- 
ate between different atomic and molecular species. Other physical processes 
(like vaporization, ionization, dissociation, excitation, chemistry, etc.) can be 
added in the future to expand the applicabilit y of the code. The Mont e Carlo 
methodology used in the code is described by Xie and Mumma ( 1996| ) for the 



case of cometary atmospheres. Our modifications are related to the physical 
conditions in the earth's atmosphere and the sputtering process. Here we de- 
scribe the basic numerical procedure, while for the detailed description we 
refer to the original paper. 

The numerical simulation follows the flight of a meteoroid through the atmo- 
sphere and calculates how many particle collisions with the air happen during 
one time step St of calculation. Each modeled particle in DSMC is a statistical 
representative of a set of particles. For N MC model particles within time step 
5t at altitude H, the number of real particles in one set is 

N set (H) = ^mn^(H)6t(H) (?) 
Nmc 



If 5t is held constant during the simulation then the distance covered by the 
meteor during one time step would eventually become much larger than the 
atmospheric mean free path Z atm . In order to avoid this statistical problem, we 
move the meteoroid one Z atm per step, thus St(H) = la,tm(H)/v m . The number 
of Monte Carlo particles used in our computations is iV MC =i500. 

A model particle is ejected in a random direction from a random location 
on the meteoroid surface. We use a spherical meteoroid of _R m =0.05 m and 
isotropic sputtering for simplicity. A more general shape and sputtering angle 
distribution could be used in the future. We also assume a constant sputtering 
velocity V^ p of 20 km/s relative to the meteoroid. This corresponds to ~10% 
of the collision energy if the masses of projectiles and sputtered particles are 
equal. We plan to introd uce a more reali stic distribution of sputtere d velocities, 
such as a Maxwellia n ( Coulsonl . 20021 ) or Sigmund-Thompson ( Thompson . 



196ilSigmun"dl . ll969h . 



The history of collisions of an ejected particle and all subsequent collisionally 
produced fast air particles is recursively traced by the code. The basic as- 
sumption is that fast particles collide only with the background atmosphere, 
but never with each other. This is a good approximation as long as the local 
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number density of ejected particles around the meteoroid is not comparable 
to the atmospheric density. This condition is satisfied in our case of sputtered 
particles. 

Since all modeled particles are treated as the same species, there is only one 
collisional cross section a. We use the Variable Hard Sphere model where the 
cross section is velocity dependent: 




We use the reference values typical for air: s=0.25, <7 re f=1.26xl0 9 m and 



values ty_ 

w re f=(4989/M)°- 5 m/s ( Bovdl . l2000j ) . where M is the molar mass (kg/mol) of 



air at the altitude of collision. The time between collisions for a particle of 
velocity v is calculated from a random number between zero and one: 



w = E At (#)- ( 9 ) 

H 



H 



This was derived from the assumption that air particles, whose velocities are 
much smaller than the velocity of modeled particles, follow the Maxwellian 
distribution. The sum over n atm (-ff) is introduced because the atmospheric 
density is changing with altitude along the particle trajectory. To describe 
this change we use altitude steps of 1 km. The time that the particle spends 
crossing an altitude layer is At(H). Since neither t co n nor H is known in 
advance, the equation is solved iteratively. 

A collision between particles requires random selection of the velocity vector 
Vt of an air particle in the Maxwellian distribution. We choose cylindrical 
coordinates (<pt,v p t,v z t) where 

v t = v pt cos <p t i + v pt sin ip t j + v zt k. (10) 
From three random numbers 9^2 and 9^3, the velocity is given as: 

Vt = 27T^! (11) 

(12) 

Erf I i^i ) = 2|9ta — 0.5) (13) 
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where Erf is the Error integral, v 2 =2RT(H)/M(H), and T(H) and M(H) are 
the temperature and molar mass of the air at altitude H . The collisional scat- 
tering of two particles is isotropic, with velocity directions chosen randomly in 
the center-of-mass frame. When the kinetic energy of a particle drops below 
some pre-defined limit (10 eV in our simulations), its trajectory and collisions 
are not traced any more. This limit cannot be below the energy of photons 
that we expect to be produced. Another reason for losing a particle is colli- 
sion with the meteoroid. We do not consider sputtering from such secondary 
collisions. 



4 Model results 



In order to visualize a meteor, we collect the total number of collisions visible 
from the ground by a TV camera at location r cam . The coordinate system is 
shown in figure 3. The meteor's direction of flight is v m = — sin Z j — cos Z k, 
where Z is the zenith angle (45° in our model). The camera's frames point at 
subsequent locations along the meteor path in order to track the meteor flight. 
When a collision happens at r co ii in space and t co ii in time (starting from the 
first frame), the camera's frame number / cam is calculated (by ignoring decimal 
places) from 

, t C o\i + l^coii _ ' cam | / ° / -i . \ 

/cam — ^ v J 



where c is the speed of light. The number of collisions in one set iV S et (equation 
7) is added to the image pixel that points at the collision. 



The first appearance of the meteor depends on the camera's sensitivity. The 
pixel values can be transformed into apparent magnitudes if we define the 
frame in which the meteor first appears. This is the frame showing the meteor 
at its beginning height i/begin- The brightest pixel in that frame is set to be 
the limiting magnitude SPtiimit and its value, iVi imit , is used for calculating the 
apparent magnitude. If the distance of the meteor from the camera in this 
first frame is Dumit then the brightness of a pixel is 



m = 2rt« m ,t - 2.5 Ioe 



-^limit \ -^collisions 

D J ~Nv 



limit 



(15) 



where iV co u is i ons is the pixel value and D is the meteor's distance. 



Figure 4 shows images of our meteor model from f our different camera l oca- 
tions. The basic morphological features described by Spurnv et al. ( 2000b[ ) are 
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successfully reproduced. The meteor first appears as a diffuse structure that 
grows in size, but then turns into a comet-like V-shape, followed by a transi- 
tion into a more compact, elongated structure. The exact shape depends on 
the relative position between the camera and the meteor path. Since we do not 
consider photon emission from long-lived atomic states, the feature missing at 
lower altitudes, below ~140 km, is the train that develops behind the meteor. 
Figure 4 also shows the number of collisions "visible" by each pixel. Collisions 
are concentrated at the center of the diffuse structure, which demonstrates 
the plausibility of the value used for \l> in equation 3. 

The efficiency of light production in high altitude meteors (described as 7 in 
equation 3) depends on the number of collisions. Collisional cascades differ in 
their total number of collisions. This number is mostly between 100 and 200, 
with ~180 collisions being the most frequent (figure 5). The total energy of 
collision between two particles follows the distribution dN = Nf(E)dE that 
can be calculated from 

AN 

= m (16) 



where AN and AE are small increments in the number of collisions and 
collisional energy, and the total number of collisions is N. The distribution 
obtained is shown in figure 6. For most of its range it has the functional 
dependence f(E) oc E~ p , with p=0.66. 

From the meteoroid's point of view, the diffuse coma does not change much 
with altitude when measured in the units of atmospheric mean free path / atm . 
Figure 7 shows the meteor viewed by a camera that moves with the meteor. 
Since collisions between the fast coma particles are not considered, the size of 
the coma is controlled exclusively by / atm - Decreasing altitude only increases 
the number of collisions per unit volume. If the assumption of negligible coma- 
coma interactions holds all the way to the end of the molecular flow regime, 
then the meteor would shrink to a few meters in size at 120 km altitude, and 
below one meter at ~100 km. The fact that the observed meteor size is larger 
at these altitudes indicates that the coma eventually becomes dense enough 
to produce internal collisions. In addition, vaporization starts to produce par- 
ticles in larger numbers than sputtering, making the coma denser than during 
sputtering. 

Meteor images based on apparent magnitudes (upper rows in figure 4) can be 
used to calculate the meteor light curve. Total meteor brightness is the sum 
of all pixels brighter than the limiting magnitude. Hence the absolute meteor 
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magnitude, denned as the brightness at 100 km distance, is 

M = 9Jl limit -2.51og(-^) 2 -2.51ogf £ 

Vl00km7 V-Pi X <6-5 m Nnmit J 



(17) 



where m pix is the pixels' apparent magnitude. Since M is a function of Amut> 
the light curve depends on the beginning height H hegin and on the camera's 
position. The dependence on i/begin is stronger, because iVi imit is a function of 
altitude. Notice that according to equation 3, iVn m i t is a constant for a given 
meteor and detector. By letting A^n m i t vary with altitude, we declare it an 
unknown that can be determined experimentally from the pixel value of the 
limiting magnitude in the first meteor frame. 

Figure 8 shows the light curves of our meteor model for two beginning heights: 
200 km and 171 km. The cameras have different distances to the starting point 
of the meteor, thus they differ in the absolute magnitude at H begin . In the first 
few frames, the meteor brightness is slightly erratic as the meteor is a few 
pixels in size and susceptible to stochastic variations in pixel values. After 
that the light curve smoothes out, but it is not linear; it shows a change 
i n slope, which is als o an observed characteristic of high altitude meteors 
( Spurnv et all l2000bh . 



5 Conclusion 



We have shown that thermalization of particles sputtered from a meteoroid 
surface by the impinging atmosphere can explain the phenomenon of high 
altitude meteors. A sputtered particle has a much larger kinetic energy than 
the surrounding atmospheric particles. It takes between 100 and 200 collisions 
between particles in the resulting collisional cascade to redistribute this initial 
energy down to the thermal level of the atmosphere. The volume filled with 
collisions surrounding the meteoroid is controlled by the atmospheric mean 
free path. These collisions are the source of the atomic excitations required 
for the light production. 

We showed that this process can create enough photons to be detected by high- 
sensitivity TV cameras. The number of sputtered particles, and consequently 
the number of photons produced, is proportional to the meteoroid size. This 
explains why the observed beginning height correlates with the maximum 
meteor brightness. Our analytical description of this correlation agrees with 
the data. We developed a Monte Carlo code specifically for the problem of 
high altitude meteors. Simulations successfully reproduced the observed size 
and morphology of these meteors. We found that their shape also depends on 
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the relative angle between the meteor path and the camera's line of sight. 



This model provides the grounds for further study of meteoroid properties 
and meteor microphysics. Features like streams and jets observed in high al- 
titude meteors are probably created by nonuniform sputtering and meteoroid 
rotation. This will be addressed in the future by Monte Carlo modeling. It 
is also possible to trace different atomic species during thermalization, which 
is a highly nonequilibrium process, and to calculate their spectral lines. Such 
numerical simulations will remove the need for "fudge" factors typically used 
to hide underlying microphysics in meteor equations. 
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APPENDIX: Maximum brightness from the single body ablation 
theory 



Meteor altitude H is described by the equation of motion: 
dH 

— = -v m cosZ. (18) 
at 

Ablation is changing the meteoroid mass as: 

^(xn atm (H)m 2 / 3 v 3 m . (19) 

Combining and integrating these two equations, while keeping the meteor 
velocity constant (a good approximation for fast meteors like Leonids) and 
assuming an exponential atmosphere, we get the solution: 

m^iH) = m l J? - constant x (20) 

cos Zi 

where m m is the initial meteoroid mass. At the end point of the trajectory 
i^end, the meteor mass becomes zero, thus 

n a tm(#end) oc m^ 3 COS Z. (21) 
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The maximum meteor light intensity is (e.g. Bronshten . 1983t Opik . 1958| ) 

/max OC m^ /3 n atm (-ffend)- (22) 



Since m^ 3 oc R m , from equation 22 and 21 it follows that the maximum 
brightness in magnitudes is 



— 7.51ogi? m — 2.5 log(cos Z) + constant. 



(23) 
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Fig. 1. Beginning heights as a function of the maxim u m brig htness for Leonid me 



teors. Diamond s are o bservations from ISpurnv et al.l (2000a) and circles are from 



Campbell etafl {2000 ). Even though these two data sets are from different obser 



vational setups, they follow a similar functional dep endence indicated by the solid 



line. Two Leonids observed bv lFuiiwara et al.l (|199a ) are shown as squares 



Fig. 2. Beginning heights as a funct ion of the reduced m aximum brightness (see 
equation 6). Only observations from Spurnv et al. 1 2000al ). marked as diamonds, 



are used, in order to have the same limiting magnitude for all data. The solid line 
is the linear regression, with the 95% confidence interyal sho wn as dashed lines. For 



comparison, two Leonids observed bv lFuiiwara et al.l (|1998l ) are shown as squares 



Fig. 3. Geometry of the modeled meteor and camera positions. The meteor's entry 
(zenith) angle is 45° in the Y-Z plane and its trajectory points toward the origin 
of the coordinate system. The cameras' coordinates are: (0, 150km, 0) for Camera 1, 
(50km,100km,0) for Camera 2, (100km, 50km, 0) for Camera 3 and (0,-100km,0) for 
Camera 4. 

Fig. 4. Modeled meteor observed by Cameras 1-4 (see figure 3), respectively, from 
top to bottom. The upper row in each panel shows the meteor brightness in apparent 
magnitudes (see equation 15) based on the beginning height of 200 km and the 
limiting magnitude of 6.5 m . The meteor's height and distance from the camera are 
indicated in each frame, together with the pixel size. The angular size of pixels is 
3.75" x 3. 75". The lower row in each panel shows the number of particle collisions 
visible at each pixel in these same camera frames. The meteoroid radius is 0.05 m. 

Fig. 5. Rate of occurrence of the total number of collisions within one collisional 
cascade, calculated from 100,000 cascades. 

Fig. 6. Distribution of collisional energies (see equation 16) during thermalization, 
calculated from 1.6xl0 7 collisions. 

Fig. 7. Modeled meteor observed by a camera that follows the meteor. The pixel 
size, indicated in each frame, is equal to the atmospheric mean free path at the 
meteoroid altitude. The asymmetry at larger altitudes comes from the variation of 
the mean free path in the coma. 

Fig. 8. Theoretical light curves of the high altitude meteor model for the beginning 
heights of 200 km and 171 km, observed from four different locations on the ground 
(figure 3). 
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Fig. 5: 




